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Abstract: Central exclusive production (CEP) processes in high-energy hadron-hadron 
collisions provide an especially clean environment in which to measure the nature and quan- 
tum numbers (in particular, the spin and parity) of new resonance states. Encouraged by 
the broad agreement between experimental measurements and theoretical predictions based 
on the Durham approach, we perform a detailed phenomenological analysis of 77 and me- 
son pair CEP final states, paying particular attention to the theoretical uncertainties in 
the predictions, including those from parton distribution functions, higher-order perturba- 
tive corrections, and non-perturbative and proton dissociation contributions. We present 
quantitative cross-section predictions for these CEP final states at the RHIC, Tevatron 
and LHC colliders. 
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1. Introduction 



There has recently been a renewal of interest in studies of central exclusive production 
(CEP) processes in high-energy proton-(anti)proton collisions (see [1-4] for reviews), both 
theoretically [5-14] and experimentally [4,15-31] (see [32,33] for reviews). The CEP of an 
object X may be written in the form 



where + signs are used to denote the presence of large rapidity gaps. An important 
advantage of these reactions is that they provide an especially clean environment in which 
to measure the nature and quantum numbers (in particular, the spin and parity) of new 
resonance states, from 'old' Standard Model mesons to Beyond the Standard Model Higgs 
bosons (see for instance [8,34-39] and references therein). One of the most interesting 
examples is the CEP of the Higgs boson, which is at the heart of the FP420 LHC upgrade 
project [40]: through the installation of dedicated forward proton detectors 420m away 
from the ATLAS and/or CMS detectors, it is hoped that detailed studies of new physics 
in high-luminosity runs of the LHC can be performed. 

As discussed in detail in [8,41,42], the CEP of, for instance, 77, heavy (c, b) quarkonia, 
new charmonium-like states or meson pairs with sufficiently large p± can serve as 'standard 
candle' processes with which we can benchmark predictions for new CEP physics at the 
LHC, as well as offering a promising way to study various aspects of QCD. In particular, 
studies of exclusive meson pair production provide a novel test of the perturbative CEP 
formalism (see e.g. [1,43]), with all its non-trivial ingredients, from the structure of the 
hard subprocess to the incorporation of various screening effects caused by absorption, 
see [44,45]. 

It is expected that a large CEP data sample will be available in the near future from 
measurements performed by STAR at RHIC [23-25], where roman pot (RP) proton de- 
tectors are ideally positioned for observing the CEP of heavy quarkonia, charmonium-like 
resonances and charged meson pairs 1 . 

As noted in [11], even without forward proton spectrometers, central diffractive pro- 
cesses of the type 



with large rapidity gaps separating the centrally produced system X from the products, 
Y and Z, of the proton (antiproton) dissociation, are still of considerable practical and 

Hn the special low-pile-up LHC runs, the TOTEM [46] and ALFA [47] detectors could study the CEP 
of hadronic systems (induced by the fusion of two soft Pomerons) using their RP proton taggers, placed at 
220m and 240m from the LHC interaction points IP5 and IP1, respectively, in combination with the central 
detectors of CMS [48] and ATLAS [49]. The expected cross sections (10-50 ^b [50,51]) should be large 
enough to allow sufficient event rates during even short runs with special LHC optics. Of particular interest 
here is the possibility of measuring azimuthal correlations between the momenta of the outgoing protons, 
which presents a way to probe the proton opacity [52]. Recently TOTEM has reported [53] encouraging 
preliminary results on central diffractive production measured with RP taggers. The possibility of measuring 
exclusive low mass pion pair production using the ALFA RP detectors during special low-luminosity LHC 
runs is also currently under discussion [31]. 



pp(p) ->• p + X + p{p) , 



(1.1) 



pp(p) ^Y + X + Z 



(1.2) 
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theoretical interest. If the incoming protons are dissociated into low masses (My t z 2 
GeV), then these reactions continue to exhibit many of the attractive properties of CEP. 
These central diffractive processes can be measured using forward rapidity gap triggers, 
with the help of simple scintillation (forward shower) counters (FSCs) [54]. Such a strategy 
was successfully implemented at the Tevatron by CDF [4], [15-19], where a rapidity gap 
trigger was used to veto on particles with pseudorapidity \rj\ < 7.4 on each side of the 
central system. The addition of FSCs to the LHC experiments could provide sufficient 
rapidity coverage, extending to the 'blind' region 6 < |r/| < 8, which would allow the 
exclusion of events with high-mass and a large fraction of events with low-mass diffractive 
dissociation in the special low pile-up LHC runs. The FSC counters proposed for CMS [54, 
55] have recently been installed in the LHC tunnel and can be operational for low pile-up 
conditions which require special LHC running modes. A proposed dedicated CMS run with 
an integrated luminosity of about lOOpb -1 of single no-pile-up interactions, with FSC (and 
Zero Degree Calorimeters (ZDC) which detect neutral particles) engaged and taking data 
in common with TOTEM is especially promising 2 . It is also worth mentioning that by 
employing FSCs a rich CEP physics programme could also be realised with the LHCb [27- 
29] experiment. The excellent particle identification and high momentum resolution of the 
LHCb detector are especially beneficial for measurements of low multiplicity final states. 
In the near future, LHCb is planning to perform measurements of exclusive meson pair 
production using the low multiplicity dihadron trigger, including Xc decays into 7r + 7r~ and 
K + K~ states [27,56]. Searches at LHCb for the baryonic (pp, AA) decay channels of the Xc 
mesons and new charmonium-like states are also now under discussion [57]. A promising 
study of low central mass CEP events is also ongoing at ALICE [26], using additional 
scintillator detectors placed on both sides of the central barrel, which allows tagging of 
double rapidity gap events. 

In 2007 CDF published a search for 77 CEP [16] at the Tevatron, with £'±(7) > 
5 GeV. Three candidate events were observed, in agreement with the expectation of [42]. 
Subsequently, to increase statistics the £±(7) threshold has been decreased to 2.5 GeV, 
and in [19] the observation of 43 77 events in [77(7)! < 1.0 with no other particles detected 
in —7.4 < 77 < 7.4 is reported, which corresponds to a cross section of <7 77 = 2.48^35 
(stat)lo5i (syst) pb. The theoretical cross section, calculated using the formalism de- 
scribed in [8,42] and implemented in the SuperCHIC MC generator [58], is 1.42 pb using 
MSTW08LO PDFs [59] and 0.35 pb using MRST99 (NLO) PDFs [60], while the P± , A<j> 
and invariant mass distributions of the 77 pair are well described by the MC. In the anal- 
ysis in [19] special attention was paid to the possible background caused by 7r°7r° CEP, 
since one or both of the photons from 7r° — > 77 decay can mimic the 'prompt' photons 
from gg — > 77 CEP. Importantly, CDF has found that the contamination caused by 
CEP is very small (< 15 events, corresponding to a ratio N(ir°7r°)/N('y"/) < 0.35, at 95% 
CL). This supports the non-trivial theoretical result of [11], calculated within a perturba- 
tive QCD framework, which predicts o-(7r°7r°)/<7(77) ps 0.01 for the CDF event selection 3 . 

2 We thank Mike Albrow for a valuable discussion about this. 

3 We recall that, due to the J z = 'selection rule' which operates for CEP, the LO gg — > tv°tt" amplitude 
vanishes for J z — incoming gluons, leading to a strong suppression in the production cross section. 
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This adds another strong argument in favour of further developing a better quantitative 
theoretical understanding of the CEP of hadron pairs (h = n, K,p, A) 

pp{p) -^p + tih + p(p) (1.3) 

in the spirit of the discussion in [11]. We recall that these channels, especially tt + tt~ and 
K + K~, are also ideally suited for the spin-parity analysis of the Xc mesons, and detailed 
studies of the CEP of such states provide a promising way to probe various aspects of 
QCD, see for example [8-11], [41]. 

Recently LHCb has reported preliminary results on the CEP of Xc mesons in the 
Xc J/ip + 7 channel, where vetoing was imposed on additional activity in the rapidity 
region 1.9 < rj < 4.9, with sensitivity to charged particles in the backwards region —4 < 
r] < —1.5 [27-29,56]. While the Xc(o,i) production data are in good agreement with the 
theoretical predictions for exclusive production [5,8], the observed Xc2 rate is somewhat 
higher. However, it is worth recalling that the observed LHCb data include some fraction 
of events with proton dissociation, as in (1.2). In [11] qualitative arguments were given 
that the protons dissociative process should favour the production of higher spin Xc(i,2) 
states, with the Xc2 yield being particularly enhanced. Recall that in the non-exclusive 
case the momentum k± transferred through the i-channel gg system (see Fig. 1) can be 
rather large, which leads to an increase of the higher spin Xc(i,2) contributions. The Xc2 
central production cross section, which is proportional to k^_, is in particular expected to 
be sizeably enhanced. However a more accurate account of the effects caused by the un- 
instrumented regions in the LHCb experiment [27-29] requires more detailed quantitative 
studies 4 . 

We note that CMS has published [62] results of a measurement at 7 TeV during the 
2010 LHC run of exclusive dimuon production mediated by two-photon fusion. This could 
be an important step towards forthcoming CEP studies with the CMS detector at the 
LHC. The observation of 77 CEP events, a search for which was recently reported by CMS 
in [63], would in particular be of much interest 5 . 

The high precision of the new CDF data on 77 CEP [19] and the expectation of new 
results from the RHIC and LHC experiments on central diffractive production encourages 
interest in a more detailed theoretical understanding of these processes, and in reducing 
the theoretical uncertainties. Apart from the dependence on the choice of the set of parton 
distribution functions (PDFs) used for the computation of the CEP cross sections, the 
main uncertainties in the theoretical predictions come from possible NLO effects 6 and 
uncertainties in the evaluation of the soft absorptive effects. In this paper we revisit these 

4 On the experimental side, the addition of FSCs on both sides of the LHCb experiment [61] would allow 
a more efficient veto on inelastic events and should greatly clarify the situation. 

5 We recall that such a measurement, and indeed measurements of other CEP processes more generally, 
will however become increasingly difficult during the current high luminosity LHC runs, as the pile-up 
rate prevents an effective use of the rapidity gap veto method. As described above, this problem may in 
principle be overcome by making use of proton taggers, which may be viable in the case of for example rnr 
CEP [31], but in the case of 77 CEP this may prove very challenging due to the quite low expected cross 
sections. 

6 At present the predictions in [5], [8-11] are based mainly on the LO formulae. 
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theoretical uncertainties, paying special attention to the structure of the NLO corrections, 
in particular for the case of 77 CEP. Based on the most recent Durham models of soft 
diffraction [50,51], we update the expectations for both the eikonal and semi-enhanced 
absorptive corrections. We also present numerical results for LHC and RHIC energies 
accounting for both the expected experimental cuts and acceptances. 

The 7T7T CEP process, mediated by Pomeron-Pomeron fusion, has been a subject of 
theoretical studies within a Regge-pole framework since the 1970s (see, for instance [64-66] 
for early references and [9, 10, 12] for more recent ones). There have also been a variety of 
experimental results on low mass meson pair CEP, in particular from the CERN-ISR, with 
the cross sections in broad agreement with the expectations from Regge phenomenology 
(see [2,4] for reviews). As discussed in [9,10], at comparatively large meson transverse 
momenta k± CEP should be dominated by the perturbative 2-gluon exchange mechanism, 
where the gg — > tttt, K + K~ subprocess can be modelled using a generalisation of the 
formalism of [67,68]. In the experimentally relevant kinematic region (M nn ~ M x , /c_i_(7r) ~ 
M x /2) in principle both (non-perturbative and perturbative) mechanisms could contribute 
to the overall rate and this issue requires careful investigation. Of special interest is the 
transition region between these two contributions which is very sensitive to the behaviour 
of the meson form factor Fm{P)- Here we perform a detailed study of the CEP of meson 
pairs paying special attention to the perturbative regime, and the transition between the 
non-perturbative and perturbative regimes. 

Finally, we note that for the sake of brevity we will not discuss the case of exclusive 
dijet production here (see for instance [69,70]). Although it is certainly an interesting 
channel through which we can test the pQCD CEP framework out to quite high Mjj, it 
is not without issues. In particular, in the dijet case it is hard to define purely exclusive 
kinematics and to exclude the contribution of events with a relatively soft third jet [70]. 
There are also additional uncertainties to consider, in particular in the value of jet's Et, 
caused by the details of the jet searching algorithm. 

This paper is organized as follows. In Section 2 we recall the main aspects of the per- 
turbative formalism for CEP. In Section 3 we consider in some detail the impact of NLO 
corrections and the choice of parton distribution functions on our CEP predictions, using 
77 CEP as a paradigm 'standard candle' process. In Section 4 we discuss the proton dis- 
sociation contribution to standard 'elastic' CEP, which is relevant for when CEP processes 
are identified by large rapidity gaps, rather than by explicit tagging of the final-state pro- 
tons. In Section 5 we expand on our previous treatment of meson-pair CEP by considering 
two other potentially important contributions, the first (perturbative) arising from the so- 
called 'symmetric' contribution and the second (non-perturbative) from double Pomeron 
exchange. In Section 6 we present and discuss our numerical predictions for meson pair 
and 77 CEP at the LHC, Tevatron and RHIC, and in Section 7 we summarize our results 
and make some concluding remarks. 

2. Central exclusive production 

The formalism used to calculate the perturbative CEP cross section is explained in detail 
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elsewhere [1,8,71,72] and so we only review the relevant aspects here. The amplitude is 
described by the diagram shown in Fig. 1, where the hard subprocess gg — > X is initiated 
by gluon-gluon fusion and the second £-channel gluon is needed to screen the colour flow 
across the rapidity gap intervals. We can write the 'bare' amplitude in the factorised 
form [1,34,35,41,42] 



7T 2 



<i 2 Q KA 

—2— ^— ■ r2fg{xi,x' 1 ,Q 2 1 ,fi 2 ;t 1 )f g (x 2 ,x' 2 ,Ql,iJ, 2 ;t2) , (2.1) 

QKQ.L - pi ± r(Q± + P2 X r 



where the / s 's in (2.1) are the skewed unintegrated gluon densities of the proton: in the 
kinematic region relevant to CEP, they are given in terms of the conventional (integrated) 
densities g(x,Q 2 ). ti is the 4-momentum transfer squared to proton % and [i is the hard 
scale of the process, taken typically to be of the order of the mass of the produced state: as 
in [35,41], we use \x = Mx/2 in what follows. The f-dependence of the / s 's is isolated in a 
proton form factor, which we take to have the phenomenological form Fj^(t) = exp(bt/2), 
with b = 4GeV -2 . M is the colour-averaged, normalised sub-amplitude for the gg — )■ X 
process 



Here a and b are colour indices, Mx is the central object mass, represents the gg — > X 
vertex and qi ± are the transverse momenta of the incoming gluons, given by 

<?l x =Q±~ Pi j. , Q2±_ = ~Q± - P2±_ > ( 2 - 3 ) 

where Q± is the momentum transferred round the gluon loop and pi ± are the transverse 
momenta of the outgoing protons. Only one transverse momentum scale is taken into 
account in (2.1) by the prescription 

Qt = min{Q_L, |(Qj_ - Pi ± )|} , 

Q 2 = min{Q ± ,|(Q ± + P2± )|} . (2.4) 



'cik 



^onh 




fg(x2,---) 

Figure 1: The perturbative mechanism for the exclusive process pp — > p + X + p, with the eikonal 
and enhanced survival factors shown symbolically. 
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The longitudinal momentum fractions carried by the gluons satisfy 



where x' is the momentum fraction of the second t-channel gluon. The differential cross 
section at X rapidity yx is then given by 

dcr lo2 \ f j2„ l r (Pi±.P2±))| 



dy^ 



(Si lh ) J d 2 p lx d 2 p 2x ^g^ g e 2 ik (p lx ,p 2 J , (2.6) 



where T is given by (2.1) and S 2 ^ is the 'eikonal' survival factor, calculated using a gen- 
eralisation of the 'two-channel eikonal' model for the elastic pp amplitude (see [52] and 
references therein for details). 

Besides the effect of eikonal screening, S ei k, there is an additional suppression caused by 
the rescatterings of the intermediate partons (inside the unintegrated gluon distribution, 
f g ). This effect is described by the so-called enhanced Reggeon diagrams and usually 
denoted as S 2 nh , see Fig. 1. The value of S 2 h depends mainly on the transverse momentum 
of the corresponding partons, that is on the argument Q 2 of f g (x, x' , Q 2 i,H 2 ;t) in (2.1), and 
depends only weakly on the p± of the outgoing protons (which formally enters only at 
NLO). While in [5,8] the S 2 nh factor was calculated using the formalism of [73], here, 
following [11], we use a newer version of the multi-Pomeron model [74] which incorporates 
the continuous dependence on Q 2 and not only three 'Pomeron components' with different 
'mean' Qi. We therefore in practice include the S eri h factor inside the integral (2.1), with 
(5 2 nh ) being its average value integrated over Q±. 

If we consider the exact limit of forward outgoing protons, pi ± = 0, then we find that 
after the Q± integration (2.2) reduces to [1] 

M OC q\ ± qi ± V tJ -> \q\{V ++ + F__) ~ SXlX2V ^2 > ( 2 - 7 ) 

where A( 12 ) are the gluon helicities in the gg rest frame. The only contributing helicity 
amplitudes are therefore those for which the gg system is in a J z = state, where the 
z-axis is defined by the direction of motion of the gluons in the gg rest frame, which, 
up to corrections of order ~ q 2 ^/M x , is aligned with the beam axis. In general, the 
outgoing protons can pick up a small p± , but large values are strongly suppressed by the 
proton form factor, and so the production of states with non-J z = quantum numbers is 
correspondingly suppressed (see [5,8,75] for examples of this in the case of X(c,b) an d V(c,b) 
CEP). In particular, we find roughly that 

|r(|j,i = 2)r.. w 



\T(J, = 0)|= (Qi> 2 ' 

which is typically of order ~ 1/50 — 1/100, depending on the central object mass, c.m.s. 
energy and choice of PDF set. As discussed in [11], this ' J z = selection rule' [34,71,76] 
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will have important consequences for the case of meson pair CEP. Finally, we note that in 
(2.7) the incoming gluon (transverse) polarizations are averaged over at the amplitude level: 
this result is in complete contrast to a standard inclusive production process where the 
amplitude squared is averaged over all gluon helicities. Eq. (2.7) can be readily generalised 
to the case of non-J z = gluons which occurs away from the forward proton limit, see in 
particular Section 4.1 (Eq. (41)) of [8], which we make use of throughout to calculate the 
MM CEP amplitudes from the corresponding gg — > MM helicity amplitudes. 



3. 77 CEP and NLO corrections 
3.1 77 CEP revisited 

As discussed in the Introduction, in [16] CDF published a search for 77 CEP with E±(j) > 
5 GeV, observing three candidate events. Subsequently, to increase statistics the E±(j) 
threshold was decreased to 2.5 GeV, and in [19] the observation of 43 candidate 77 events 
in 177(7) I < 1.0 with no other particles detected in —7.4 < 77 < 7.4 was reported, which 
corresponds to a cross section of <r 77 = 2.48t(X35(stat)t(X5i(syst) pb, and with a 7r°7r° 
contamination consistent with zero. The theoretical cross section, calculated using the 
formalism described in [8,42] and implemented in the SuperCHIC MC generator [58], is 
1.42 pb using MSTW08LO PDFs [59] and 0.35 pb using MRST99 (NLO) PDFs [60]. Evi- 
dently the prediction using the LO PDF set is consistent with the result within theoretical 
uncertainties, although both predictions lie below the observed cross section, particularly 
so in the case of the MRST99 PDF choice. 

It is natural to ask why the 77 CEP cross section predictions in [8] are somewhat lower 
than the data. In fact, there are reasons why we may expect this to be the case. Most 
importantly 7 , we recall that the prediction of [8] includes only the LO perturbative QCD 
contribution to the 77 CEP process. In general, we may reasonably expect a numerically 
large NLO K-factor correction to the gg — > X subprocess: for example, the higher order 
corrections to Standard Model Higgs boson production via gg — >■ H (see [77, 78] and refer- 
ences therein) and P-wave quarkonia decay X ^ 99 [79] are known to be quite large 8 . In 
the case of inclusive Higgs boson production via gg fusion, the K-factor is approximately 
given by 



K » 1 + 



a s (Mfj 



TT 



2 11 



1.5 , (3.1) 



for a light Higgs (Mh ~ 125 GeV). We can see in particular that the NLO cross section is 
numerically enhanced by a factor of tt 2 , the origin of which can be traced back to a Sudakov- 
like double logarithm, as log 2 (— \q 2 \/Mfj), when the imaginary part of the logarithm is 
squared. Such an enhancement is also observed to occur in the NLO correction to the 



Experimentally we may also expect the observed 77 cross section to be enhanced by the small fraction 
of 77 events seen by CDF which are not truly exclusive, but rather due to double diffractive production 
where one or both of the proton and antiproton dissociates, but where the proton dissociation products are 
not seen within the CDF acceptance; we estimate such a fraction to be 10% or lower. 

8 In the case of Xc,b CEP we assume the same K-factor for the x 99 an d gg — >• X cases, see [5,8], 
which is only true to a certain degree of approximation. 
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X(c,6)0 — > 55 decay. We may reasonably expect similarly enhanced NLO corrections to 
Higgs and X(c,b)o CEP, and in [1,5,41] such an assumption is made to estimate the relevant 
NLO cross sections, although we note that the K-factor in (3.1) is for the spin-averaged 
case, whereas in CEP only the J z = helicity amplitudes are important. It is also worth 
mentioning that an additional positive correction (of the order of about +20%) could come 
from a consistent account of the self-energy insertions in the propagator of the screening 
gluon, in the spirit of the discussion in [80]. 

In the case of 77 CEP, the NLO corrections to the g(\i)g(\2) — > 7^3)7^4) subprocess 
helicity amplitudes are calculated in [81]. While it is quite difficult to trace the origin of all 
the 7r 2 terms in the reasonably complicated analytic expressions for the helicity amplitudes 
in [81], our understanding is that the 77 NLO virtual corrections may be similarly enhanced, 
e.g. by Sudakov-like double logarithmic terms. Indeed, including the finite part of these 
NLO corrections [81] to the gg — > 77 subprocess, we find that the K-factor for the 77 CEP 
cross section is quite large, given by approximately K ~ 1.6 for experimentally relevant 
M 77 values. However, as discussed in the following section, this does not constitute the 
full NLO K-factor, and therefore can only be taken as a guide. Nonetheless, this result is 
suggestive that the LO estimate for the 77 CEP cross section may receive sizeable positive 
higher order corrections. 

This effect could easily give a factor of ~ 2 increase in the theoretical 77 CEP cross 
section, when compared to the calculation in [8], thus alleviating some of the tension that 
exists between the theory and the latest CDF measurement. Nonetheless, in the case of 
the MRST99 PDF choice, the theoretical prediction is almost an order of magnitude below 
the experimental value, and so it appears that, even taking into account the other sizeable 
theoretical uncertainties as well as the issues discussed above, the 77 CEP cross section 
calculated with this PDF choice is too small. This may indicate that in fact the 'true' 
gluon PDF in the low-x, low-Q 2 region relevant to CEP is better described by the LO 
MSTW08 fit, which gives a higher density of gluons due to the steeper x dependence than 
the relatively flat MRST99 PDF. Of course, given the existing theoretical uncertainties 
in the CEP calculation, it is difficult to make a very strong statement about this, and 
certainly other measurements, in particular of 77 (and other processes) at the LHC, would 
be very useful in shedding light on this issue. 

We note that the proton's PDFs are physical quantities which are defined as the 
expectation values of the corresponding QCD operators over the proton state. There are 
no indices, such as LO, NLO, MS,... , in their definition. However, these PDFs are not 
measured directly, rather they are determined by optimising the description of some hard 
cross section by the convolution of the PDF with the corresponding coefficient function 
(i.e. the hard matrix element |.M| 2 ). At this stage the resulting PDF starts to depend 
on the perturbative approximation (LO, NLO,...) and the choice of factorisation and 
renormalisation schemes used for the calculation of \M\ 2 . For example, for the Drell-Yan 
process the K-factor K = | ^A^ 1 ^ 1 " 1 2 / 1 7W LO | 2 ~ 1.5 — 2 is rather large, which will lead 
to a larger PDF using the LO than when using the NLO \M\. However we cannot 
claim that different processes have the same K-factors. For example, in 77 CEP the NLO 
correction is expected to be numerically of the same order (for the hard |.M| 2 ) as that in 
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the DY case but in CEP the value of \Ai\ 2 is multiplied by (xg) . 

In practice, global fits of the gluon density, g(x, Q 2 ), at low x are based almost entirely 
on measurements of F 2 {x, Q 2 ), in particular through the effect that g(x,Q 2 ) has on the 
DGLAP evolution of the quark PDF and therefore on dF 2 (x, Q 2 )/dlnQ 2 . The extracted 
gluon density therefore depends on the approximations used to calculate this evolution, 
but at sufficiently low scales the higher order and power (including the absorptive) cor- 
rections to the pure DGLAP evolution are certainly not negligible. In most global PDF 
analyses (performed within the linear, pure twist 2 perturbative DGLAP approach) cuts 
are imposed on the data (for example in Q 2 and/or W 2 ) in order to minimise the effect 
of these corrections, but if there is some residual contamination then this will simply be 
absorbed into the input gluon distribution. In addition, the input g{x,Qfa) is adjusted at 
low x in order to obtain a satisfactory fit to data over a wide Q 2 and x range, over most of 
which the higher-twist and higher-order corrections are expected to be small. Although 
this approach may therefore provide a satisfactory description of the gluon PDF at large 
Q 2 , it may give unreliable results for very low Q 2 and x. Indeed in the case of modern NLO 
MS PDF sets, the fits tend to prefer a negative gluon, for example g(x, Qq = 1 GeV 2 ) < 
for x < 1CP 2 in the case of MSTW08 NLO. On the other hand there are some indications, 
for instance from J/tp diffractive photoproduction data at HERA (which is not currently 
used in global PDF analyses), that the 'true' g(x,Q 2 ) at low x and Q 2 ~ 2 — 3 GeV 2 is 
larger than the current NLO PDFs [82]. 

In contrast, the LO PDFs, which generally have a g(x, Qq) which is positive at small x, 
do not give a particularly good description of the HERA F 2 (x, Q 2 ) data, and may in general 
be too large at low x and Q 2 in order to compensate for the non-negligible higher-order 
terms which are absent in the LO DGLAP evolution of the quark PDF. 9 

It is therefore not unreasonable to expect that the 'true' gluon PDF at low x and 
Q 2 , relevant to the relatively low-mass CEP processes that we are interested in here, 
lies somewhere between the lower and upper bounds set by the NLO and LO PDF sets, 
respectively. For this reason, we use the MSTW08 LO and MRST99 NLO sets in our 
study, arguing that they span a realistic range of small x gluon distributions (we show this 
explicitly in Section 6.2). Of the available NLO PDFs, we have chosen to use the older 
MRST99 set as our benchmark as this has a more benign small-x form (in particular, it 
is somewhat flat at low x, Q 2 ), and will therefore not produce numerical instabilities that 
could occur with more modern NLO fits, in which the gluon can be negative in this region 
(at least in the MS factorisation scheme generally used in the fits). 

3.2 Remarks on NLO corrections to CEP processes 

Considering the inclusive Higgs production particular example, the K-factor given 

in (3.1), due to the finite part of the virtual corrections to the gg — >■ H process, does not 
give the full NLO correction. For the fully inclusive Higgs production cross section we must 
also include the finite contribution from QCD radiation gg — > Hg (as well as qq —> Hg and 
qg — > Hq) to get the correct total NLO correction. The net effect of all these contributions 

9 In particular, at LO there is no 1/z singularity in the P qg (z) and P qq (z) splitting function, which is 
present at NLO and higher orders. 
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Figure 2: Representative NLO virtual correction diagrams to the 77 CEP process. 



can in principle be quite different from the value given in (3.1), although in fact the fully 
inclusive cross section does receive a K-factor of a similar (large) size to this. 

In the exclusive case the set of NLO diagrams that we must calculate is different. 
There are still the standard virtual NLO corrections to the 'hard' gg — > X matrix element 
(X = H,jj...), which can be calculated in the usual way or taken from known NLO 
results, such as (3.1) in the case of Higgs production. However, in CEP the contributions 
due to real emission (which cannot occur) are no longer present, and thus there is no 
explicit cancellation of the infrared divergences present in the virtual loop amplitudes. 
Here, these are cancelled by a corresponding divergence hidden in the unintegrated parton 
distributions, f g (see (2.1)). In particular, such a cancellation determines the lower k± 
limit in the Sudakov factor integral, 

T 9 (Qi,M 2 )=exp(- r 2 ^^ r A UP 99 ^) + E P ^)]dz) ■ (3.2) 



This occurs because some of the NLO loop diagrams in the hard matrix element are identical 
to those contributing to the Sudakov factor, included in the / g 's. In particular, in the 
Feynman gauge one of the first, order a s , diagrams contributing to the Sudakov factor is 
due to gluon exchange between the two incoming partons in the hard subprocess, shown 
in Fig. 2 (right) for the case of 77 CEP. As the NLO matrix element contains exactly the 
same diagram, to avoid double counting we must exclude this contribution from the NLO 
matrix element. This is just the infrared divergent contribution caused by the 'absence' of 
soft gluon emission, and thus in the NLO matrix element we should only include the region 
where the scale is larger than the argument Q of f g (Q, ■ ■ ■ ), while the 'infrared sensitive' 
low-A;_|_ part is included in the unintegrated parton densities, where there is a cancellation 
between diagrams in which the low k± gluon is absorbed by the active (right) i-channel 
gluon and the soft (left) screening gluon in Fig. 1; the long-wavelength low k± < Q± gluon 
is only sensitive to the whole colour charge of the t-channel colourless two-gluon state. As 
was shown in [43], for large subenergies squared s = M\ 3> i.e. in the BFKL limit, 
such a cancellation may be accounted for by taking the lower limit k± > Q± in the k± 
integral for the Sudakov form factor. 

In the exclusive case, we must explicitly calculate these additional diagrams, the di- 
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vergent parts of which are included in the f g 's (in particular, the 'Sudakov' loop diagram 
described above, and diagrams where additional gluons couple to the quark lines inside 
the hard matrix element, see for example Fig. 2 (left)), as well as those in which the new 
(additional) gluon couples to the 'soft' screening gluon. These will all give finite contribu- 
tions to the CEP NLO K-factor which will in general be different in the Higgs case from 
those contributing to (3.1), and these will also depend on the object, X, being produced. 
While the single and double logarithmically enhanced contributions from the higher order 
corrections to the CEP amplitude are already included in the Sudakov factor, the finite con- 
tributions coming from a full NLO calculation are not, and such a calculation has not yet 
been performed (although a subset of the contributing diagrams, in particular those which 
are logarithmically enhanced and therefore give a non-zero contribution to the Sudakov 
factor, were calculated in the Higgs case in [7]). 

Finally, a few more comments about the NLO corrections to the CEP process are in 
order. First, we note that in an NLO calculation the unintegrated parton densities, f g , 
used in (2.1) should be calculated to NLO accuracy (as described in [83]). Secondly, as 
mentioned above, there are contributions to the NLO cross section from diagrams in which 
the additional gluon is emitted from the 'soft' screening gluon but absorbed somewhere in 
the 'hard' matrix element: a representative diagram in the case of 77 CEP is shown in 
Fig. 2 (left). However, in general, when the hard scale fi ~ Mx is large, such a contribution 
should be treated as a power correction, which can be seen as follows (see also the discussion 
in [84]). We will consider two extreme cases: when the additional gluon emitted from the 
screening (left) gluon is rather soft and has transverse momentum q± ~ Q± <C Mx and 
when the gluon has a large q± ~ p. In the former case, the phase volume in the / dq\ 
integral (~ Q 2 ± ) cannot compensate the large virtuality of the additional propagator in the 
'hard' matrix element, which absorbs the gluon, and the expected contribution will be of 
the order of Q\j [i 2 . In the latter case, the new gluon has a large q± ~ fi, which would 
mean that the scale and the transverse momenta of the unintegrated partons is also large, 
with the value of q± ~ p compensated by some other parton in f g . On the other hand the 
integral (2.1) over Q± is superconvergent, i.e. it has the form 



Normally this has a saddle point at Q± = Q sp <C /i. Thus the contribution from the 
Q± ~ fi domain will be suppressed by the power factor Q 2 p / fi 2 . 

There is also some correction caused by the fact that in the 'hard' matrix element 
the incoming gluons are off-mass-shell. This also leads to a power suppressed correction 
of order Q\j since the virtuality of the incoming gluon is given by Q\ ~ Q 2 p <C p? . 
However there is some contribution from the region Q± ~ fi which, strictly speaking, can 
not be considered separately from the NLO calculation. This kinematical effect should be 
considered together with all other NLO diagrams, especially the diagrams which describe 
the creation of a highly off-mass-shell gluon from a previously low-virtuality parton, which 
may be treated as on-mass-shell. We should recall in particular that the matrix element 




(3.3) 
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with off-shell gluons is not gauge invariant and may contain some additional non-physical 
singularities which are cancelled in the sum of the diagrams which forms the gauge invariant 
group (or subgroup) of diagrams. 

In processes with large Mx 3> Q sp , such as Higgs boson CEP, these power corrections 
are numerically small and may be neglected. However at relatively low scales, as is the 
case for Xc production or 77 production when E±( n f) is not too large (as in the CDF 
event selection), these corrections may give a non-negligible contribution, although without 
performing an explicit NLO calculation it is difficult to quantify this statement. We will 
return to this issue, in particular that of the initial gluon off-shellness, in the case of meson 
pair CEP in Section 5.1. 

4. On the role of proton dissociation 

In current LHC measurements of central production only some restricted large rapidity 
gap (LRG) conditions can be imposed, leaving certain regions in the forward (backward) 
directions uninstrumented. As discussed in the Introduction (see also [11]), without forward 
proton spectrometers there is always a possibility that an 'exclusive' event does not come 
from the purely elastic process (1.1) but rather from central diffractive production (1.2), 
with sufficiently large rapidity gaps separating the system X from the products, Y and 
Z, of the proton dissociation 10 . We shall use the nomenclature Pseudo-CEP (PCEP) to 
denote such central production with the dissociation of (one or two) incoming protons. 

To evaluate the possible admixture due to proton dissociation we can use a formalism 
similar to that described in Section 2 for the CEP case, where the outgoing protons remain 
intact. To achieve this we must replace the unintegrated gluon density f g (xi, x' i: Q 2 , /U 2 ; t) 
in (2.1) by another function f g (xi,x' i ,Q 2 , fj, 2 ;t; My), which accounts for the proton dis- 
sociation into a state with mass My. The problem is that the form of this function is 
practically unknown, and so the best we can do is to make some plausible assumptions 
about its behaviour. 

For the case of low mass dissociation (p — > N*), it is reasonable to assume the same 
x, Q 2 , \i and t behaviour as for the 'elastic', p — > p, process. That is, to incorporate 
low mass dissociation it is sufficient to multiply the original (purely exclusive) result by 
some factor 1 + c, where c is the probability of proton dissociation into the relatively low 
mass state. These probabilities were measured at lower (fixed target and CERN-ISR) 
energies (see [85], and [86,87] for reviews), and can be extrapolated to LHC energies by 
accounting for stronger absorptive effects, using the models for soft diffraction discussed 
in [73,74]. This gives c ~ 0.15 for the overall probability of proton dissociation to states 
with mass My < 2 GeV, in agreement with the analysis of [85]. A similar value of c ~ 0.2 
was found at HERA, by comparing the size of the diffractive DIS cross section measured 
using the leading proton spectrometer with that found by requiring a LRG in the forward 
direction [88-90]. Thus to account for the possibility of low mass dissociation of both 

10 Note that as shown in [54,55], in the low pile-up runs ZDC detectors would allow a reduction in the 
contribution from diffraction dissociation by at least a factor of 2. 
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protons we have to increase the CEP prediction by a factor (1 + c) 2 ~ 1.4, where we take 
c = 0.2 for definiteness. 

The situation with high mass dissociation is more complicated. Usually, dissociation 
into high mass states is described in terms of triple-Pomeron diagrams. This contribution 
is driven by the triple-Pomeron vertex, g^p (see for example [86] for a review). For fixed 
momentum transferred through the Pomeron, t = —p\, that is through the i -channel state 
formed by the gluons x and x' in Fig. 1, the ratio of the high mass dissociation to the 
'elastic' contribution can be written as 



where the integral runs over the My interval available in a given experiment. The value of 
the proton-Pomeron vertex, gN(t), is known from elastic pp scattering, while the triple- 
Pomeron vertex, g^p, can be extracted from the description of lower energy (CERN-ISR, 
Tevatron) data in the triple-Regge region, see for example [91]. After accounting for the 
eikonal-like absorptive corrections, the recent fit in [92] gives <?3p(0) = 0.2^(0). However 
here we face a problem, as we will now discuss. Absorptive effects strongly depend on 
the shape of the triple-Pomeron amplitude in impact parameter, bt, space. At lower bt the 
proton optical density is large, and absorption is much stronger. On the other hand, we only 
poorly know the £— slope, &3P (that is the bt size), of the g%p oc exp(b^pt) vertex. Fits to the 
lower energy triple-Regge data [91,92] indicate that the slope b^p is rather small, and may 
even be consistent with zero, with b^p < 2 GeV~ 2 . Thus the value of the corresponding 
survival factor S 2 is uncertain. Moreover, when we do not detect the outgoing proton (or 
the My state) we have to integrate over the squared momentum transfer t. In the 'elastic' 
CEP case the integral is limited by the proton form factor, and the average p± is small. 
In 'soft' single dissociation (corresponding to the triple-Pomeron diagram without any 
'hard' subprocess), the squared momentum t goes through the LRG and is again limited 
by the proton form factor. This is not the case for proton dissociation in PCEP where 
the momentum p± goes across the loop between the 'hard' subprocess matrix element 
|.M| 2 , which may be viewed as a 'point -like' blob of a very small size in bt space, and the 
g%P vertex. Now, according to (4.1), in the case of high mass dissociation we can arrive 
at very large values of p\ ~ l/bsp- This will lead to an unacceptably large probability 
for dissociation, as for large p± we cannot justify the factorization f g (x, .../J, 2 ;t, M Y ) = 
G(t)f g (x, ...ix 2 ]M Y ) of the t dependencies. Simultaneously, large values of p± will allow an 
increasing violation of the J z = selection rule which operates for pure CEP; the admixture 
of the \J Z \ =2 contribution (2.8) will increase like ~ ip]) 2 . This may be crucially important 
f° r Xc2 or 7T7T production where the J z = CEP amplitudes are strongly suppressed (the 
J z = amplitude is zero for \c2 CEP in the non-relativistic quarkonium approximation, 
and zero at LO for vrvr CEP) 11 . 

11 As was mentioned in [11], an admixture due to a proton dissociation contribution could explain, at 
least in part, why the observed LHCb [27] Xc2 PCEP cross section is higher than the CEP prediction, 
although cuts on the central system p± < 0.9 GeV are imposed to suppress this contribution. With 

the increased statistics that should hopefully be available from future data, a closer study of the dependence 




(4.1) 
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Single 


Low M Y ,z (< 2.5 GeV) 


High M Y ,z (> 2.5 GeV) 




S' 2 


0.86 ± 0.03 


0.81 ±0.03 


Double 


(Low My, Low Mz) 


(Low My, High M z ) 


(High My, High M z ) 


S' 2 


0.3 - 0.45 


0.2 - 0.28 


0.08 - 0.16 



Table 1: Estimates for the survival factor (due to strong interactions between the colliding protons) 
for single (pp — > p(Y) + 77 + Z(p)) and double (pp — >• Y + 77 + Z) diffractive two-photon 
production of a lcpton pair, at y/s = 7 TcV, for £a(lepton) > 5 GeV and for the veto range of 
1 77 1 < 5.2. The results are presented for different regions of the invariant mass, My,Mz, of the 
proton dissociation products. 



Taking b^p = 1 GeV~ 2 we then integrate (4.1) up to the maximum My allowed by the 
uninstrumented Ay for the relevant detector, and we evaluate the admixture of high mass 
dissociation in PCEP events selected by a rapidity gap veto to be about C ~ 30 — 40% 
for the CMS/ATLAS experiment and up to C — 50% in the LHCb case 12 . In other words, 
the enhancement factor by which we have to multiple the CEP cross section in order to 
estimate the PCEP rate is ~ (1 + c + C/2) 2 , where the factor of T/2' accounts for the fact 
that C is the probability that either proton dissociates. We should recall, however, that 
there is a large uncertainty in this estimate, as discussed above 13 , and that this high mass 
dissociation will generate a rather large transverse momenta p± of the central system and, 
therefore, a larger admixture of the \J Z \ = 2 contribution. 

To make the interpretation of CEP experimental data less ambiguous it would be 
helpful to suppress proton dissociation. This can be done by imposing additional veto 
conditions: no signal in the ZDC and/or, as discussed in the Introduction, FSCs, which 
cover a larger rapidity interval [54,55]. Another possibility is to select events in which 
the centrally produced system has a low transverse momentum. In particular, we can 
impose a cut on the coplanarity of the 77 or tttt pair in the transverse plane. In the loop 
formed by two Pomerons between the triple-Pomeron vertex and the matrix element of a 
'hard' subprocess, \M\ 2 , the large momentum p± goes through the \M\ 2 blob. That is, 
the centrally produced system obtains a large transverse momentum which violates the 
coplanarity of the event. By selecting events with, say, \A<f> — ir\ < 0.1 we will introduce a 
cut of P± < 0.1£_l 14 . 

Finally, let us make a few comments concerning the evaluation of the proton disso- 
ciation contribution in the case when the central system (e.g. a lepton pair) is produced 
via the photon-photon fusion mechanism. As a rule (see e.g. [62,63]), an admixture of the 
processes with proton dissociation is evaluated using the LPAIR event generator [93,94]. 
However, it is based on purely photon exchanges, and absorption effects caused by the 

of the relative Xc cross sections on this p± cut would be possible, shedding further light on this issue. 

12 In the limit that the 'hard' subprocess has characteristic E± S> 1 GeV (63P * E\ ^> 1) this estimate 
is process independent. However, for e.g. XcO production the gg — > XcO vertex may effectively enlarge the 
total i-slope in the loop somewhat, leading to a smaller probability of dissociation. 

13 Further modification of the admixture of dissociation will arise after accounting for detector effects. 

14 Recall that the A(f> distribution of CEP diphoton events observed by CDF [19], where the LRG veto 
is quite good, thus providing a very small admixture of proton dissociation events, has a narrow peak at 
A(j> = 7r, consistent with the E± of the recoil proton. 
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strong interactions between the protons are not accounted for 15 . This could therefore lead 
to some overestimation of the admixture caused by proton dissociation, while not affecting 
the purely exclusive case when the outgoing protons remain intact, where the size of the 
absorbtive effects are expected to be very small [95]. Using a technique similar to that 
discussed above for the evaluation of the proton dissociative contributions, for example, in 
the topical case of dilepton production we can evaluate the survival factors S 2 (by which 
the LPAIR predictions should be multiplied). For illustration, in Table 1 we present esti- 
mates for the case of 7 TeV proton-proton collisions with E± (lepton) > 5 GeV, and for 
the veto range \r]\ < 5.2. 



5. Meson pair CEP 

5.1 Symmetric contribution 




Figure 3: (a) Representative perturbative diagram for 'skewed' meson pair CEP. (b) Represen- 
tative perturbative diagram for 'symmetric' meson pair CEP. In both cases, quarks with momenta 
and (hih) are collinear and form colourless meson states. 



As discussed in [11], apart from the standard CEP diagram shown in Fig. 3(a), there 
is a second perturbative mechanism for producing meson pairs at high k±. This is shown 
in Fig. 3(b): the second t-channel gluon now couples directly to a quark line, while the 
collinear qq pairs can be in both flavour-nonsinglet and flavour-singlet combinations. Fol- 
lowing [11], we label Figs. 3(a) and (b) as 'skewed' and 'symmetric' CEP, respectively. 
In Fig. 3(a) the second t-channel gluon is much softer than the fusing gluons, whereas in 
Fig. 3(b) both t-channel gluons participate symmetrically in the hard subprocess. The sym- 
metric diagram represents the perturbative tail to the non-perturbative double Pomeron 
exchange mechanism, resolving the two-gluon structure of the exchanged Pomerons. 

For the calculation of the skewed contribution we use the 'hard exclusive' formalism [67] 
(see also [96]), modified to the case of gg —> MM. The amplitude can be written in the 
form ^ 

A*AiAa(M)= / dzdz'4>M(z)4^(z')T XlX2 {z,z'-rs,e) , (5.1) 
J 

15 We thank Sergey Baranov and Wenbo Li for useful discussions concerning the LPAIR MC. 



- 15 - 



where § is the MM invariant mass squared, z, z' are the meson momentum fractions carried 
by the quarks, Ai, A2 are the gluon helicities and 8 is the scattering angle in the gg c.m.s. 
frame. T\ x \ 2 is the hard scattering amplitude for the parton level process gg — > qqqq, 
where each (massless) qq pair is collinear and has the appropriate colour, spin, and flavour 
content projected out to form the corresponding meson. In the meson rest frame, the 
relative motion of the quark and antiquark is small: thus for a meson produced with large 
momentum, \k\, we can neglect the transverse component of the quark momentum, q, with 
respect to k, and simply write q = zk in the calculation of T\ 1 \ 2 - <f>(z) is the meson 
wavefunction, representing the probability amplitude for finding a valence parton in the 
meson carrying a longitudinal momentum fraction z of the meson's momentum, integrated 
up to the scale Q over the quark transverse momentum q~± (with respect to pion momentum 
k). We then use the formalism outlined in Section 2 to calculate the MM CEP cross section 
from these subprocess helicity amplitudes, see in particular Section 4.1 (Eq. (41)) of [8]. 



In [11] it was shown that the symmetric amplitude represents a power correction to 
the standard skewed CEP amplitude T, in particular with 

^ ~ ^ , (5-2) 

-'skew. rv_\_ 

where k± is the meson transverse momentum, and is the average gluon transverse 

momentum squared in Fig. 3(a): at high meson k±, the symmetric contribution will be 
subleading. However, at lower values of k± ~ Q± this will not necessarily be the case, and 
care must be taken to perform an explicit numerical evaluation of the relative contributions 
of the skewed and symmetric amplitudes. 



The symmetric CEP amplitude, considering for simplicity the case of forward outgoing 
protons (j>i ± = 0), can be written in the form 



/cPq±± 
4 X ^sym. f g (x 1 ,x l ,qj ± ,fi 2 ;t l )fg(x2,X2,qi ± ,fJ 2 ;t2) , (5.3) 

where the notation follows from Fig. 3(b), and this can readily be generalised to the case 
of non-zero proton p±. The subprocess amplitude M sy m. is given by 

Xsym. = ^^^^S b V 1± q^ULvS (5-4) 
4 

= jj^r^W • (5-5) 

For example, considering the specific case of scalar non-singlet meson production, as shown 
in Fig. 3(b), if we identify 

l\ = z'ks Z2 = (1 — z)ki I3 = zk^ U = (1 — z')ks , (5.6) 
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then we have 



y sym . = i6vr y dzdz n z )n z ) • ( 5 - 7 ) 

The other diagrams corresponding to interchanging the outgoing quark legs l\ I2 and 
'3 ^ ^4) as well as the s-channel diagrams containing 3-gluon vertices, should also be 
included, and the resulting expression can then be combined with (5.3) to give an explicit 
evaluation of the full symmetric CEP amplitude. 

Since we do not know the form of the generalized PDFs, f g , in the kinematic regime 
relevant to the symmetric CEP contribution the best we can do, following [11], is to put 
an upper limit on the cross section, based on the Schwarz inequality [97], taking 

f g (x,x',Q 2 ,...) < I((/ 9 (x,-x,Q 2 ,...) + / 9 (x',-x / ,Q 2 ,...)) , (5.8) 

where the f g (x,—x,Q 2 ,...) are the diagonal gluon distribution functions which can be 
extracted from DIS data. We note that in the symmetric amplitude we must be careful 
to forbid perturbative emission from both f-channel exchanges, which can both be hard, 
although at a scale that is somewhat lower than = Mx/2. In general we can therefore 
take 

\J Tg(Q 2 ,^i)T g (Q 2 ,fil)(xg(x, Q 2 ) + x'g(x\ Q 2 )) , 

(5.9) 

where in the case for example of Fig. 3(b), with the momenta defined as in (5.6), we have 
/if = z'(l — z)M x /4 and /j, 2 = z{\ — z')M x /4. However for simplicity and to give a rough 
estimate, we set fi± = ji2 = corresponding to (z,z') = 1/2. It should be noted, 

however, that in the region of high and low z, z' , the 'hard' scale can be arbitrarily small, 
and so the symmetric CEP amplitude may become sensitive to non-perturbative physics. 
On the other hand, the size of these corrections should be small because the wavefunctions 
4>(z), 4>{z') suppress contributions from the the endpoint regions z, z' = 0, 1. 

We also note that the expressions for the symmetric amplitude (5.4) and (5.7) include 
a non-negligible contribution from the region where the incoming gluons are far off-mass- 
shell Q 2 ~ M\ in the amplitude evaluation. For example, for the diagram shown in 
Fig. 3(b), we have 

Q2± = W ~ (1 " z )]k± ~ qi± , (5-10) 

and therefore Q\ ~ q\± ~ k 2 ^, which is of order the hard scale of the subprocess. In general, 
the incoming gluon off-shellness, or equivalently q± with respect to the beam direction, is 
generated by higher order diagrams which describe the creation of a highly virtual gluon 
from a low-virtuality incoming parton, see the discussion below (3.3). In the case where 
the dominant contribution comes from the strongly ordered region (q 2 ^_ <C M^) as in, for 
example, Higgs CEP, we can factor this initial-state gluon virtuality, generated by these 
higher order diagrams, into the q± -dependent / g 's. As we have Mx 3> Q±, corrections 
due to the off-shellness of the incoming gluons in the gg — >■ X subprocess will be small, 



f g (x, x' , Q 2 , /j, 2 , H2) 



1 d_ 

2 d log Q 2 
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and can be neglected. However, in the situation described above, this is no longer the case 
and we should in general consider a full NLO calculation where we explicitly include in the 
hard process diagrams which generate this high virtuality in the incoming gluons. Such 
a calculation may include crucial cancellations between, and contributions from, different 
NLO diagrams which would be missed simply by taking the LO expression as in (5.7), with 
the / fl 's evaluated at high gluon scales Q 2 ~ M x . We can, therefore, only consider (5.7) 
as an estimate of the expected rate: in particular the true all-orders cross section may in 
principle be somewhat larger. 

Such an issue also arises, for example, in the case of skewed tttt CEP, which we recall 
is strongly suppressed by the J z = selection rule due to the fact that the LO on-shell 
gg —7- tttt amplitudes vanish for J z = incoming gluons. Here the only non-zero contribu- 
tion arises from the case when the gluons are off-shell and therefore allow a non-J z = 
component to the cross section: for lower values of Mx, a significant contribution can come 
from the region where the gluon off-shellness is of the same order as the hard scale. The 
factorised form (2.1) we use to calculate the tttt cross section, where we ignore the gluon 
off-shellness in the gg — > tttt amplitudes, can therefore only be considered as an estimate of 
the full cross section in the lower Mx region: to give a precise evaluation, we would have 
to calculate the full NLO amplitude, where the off-shellness of the incoming gluons is gen- 
erated by additional virtual gluons. Naively taking (2.1) but with the gluon off-shellness 
explicitly included in the gg — > tttt amplitude will tend to overestimate the cross section, 
as it will miss the cancellations between different NLO diagrams that we may in general 
expect to occur 16 , as is observed in the LO gg — > tttt case, leading to the vanishing of the 
J z = amplitude, and the presence of a radiation zero in the \J Z \ = 2 amplitude. 

Finally, using (5.3) and (5.9) we can, for example, calculate the symmetric contribution 
to 7r°7r° CEP. This is plotted along with the standard skewed 7r°7r° CEP cross section in 
Fig. 4, as a function of the cut on pion E±. Also shown is the differential cross section 
da/dM^. We can see that, even at lower values of the pion E±, the symmetric contribution 
is expected to be very small, and well within other theoretical uncertainties on the skewed 
cross section prediction. As we would expect, recalling that the symmetric contribution 
represents a power correction to the skewed contribution, as E± increases this suppression 
becomes stronger. We have therefore explicitly shown that the skewed 7r°7r° contribution is 
negligible in the entire kinematic region, and this will be equally true for the CEP of other 
flavour non-singlet states (tt + tt~, K + K~ ...). Moreover, as (5.9) gives an upper bound on 
the / 5 's, the symmetric cross section may in general be even smaller than our numerical 
estimate suggests (although it is also true that a proper treatment of the gluon off-shellness 
may increase the cross section). Recalling that the tttt skewed CEP cross section is already 
strongly reduced by the J z = selection rule, this relative suppression of the symmetric 
CEP contribution will be even stronger for other meson pair production processes (rfrf, 
777/ '...) which have non-zero J z = production cross sections at LO. We can therefore 
safely neglect this symmetric contribution when it comes to producing numerical estimates 
for the expected meson pair CEP cross sections in Section 6.1. 

16 This raises doubts regarding the consistency of the strategy adopted in [12, 13], based on using explicit 
results for the amplitudes with off-shell gluons. 
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o(E± > E cut ) [pb], MSTW08LO jgj [pb/GeV], E ± > 2.5 GeV, |r;»| < 1.0, MSTW08LO 




Scut [GeV] M x [GeV] 

Figure 4: Skewed and symmetric contributions to the 7r°7r° CEP cross section for different cuts 
E cut on the tt° transverse energy, E±, and the differential cross section da/dM n7r for i?j_(7r°) > 2.5 
GeV, at -y/s = 7 TeV. In both cases a pseudorapidity cut |ry(7T°)| < I is applied. 



5.2 Non perturbative process 



Pi 




Figure 5: Representative diagram for the non-perturbative meson pair (M3, M4) CEP mechanism, 
where M* is an intermediate off-shell meson of type M. Eikonal and (an example of) enhanced 
screening effects are indicated by the shaded areas. 

The formalism of Section 2 is only valid for relatively large values of the central system 
mass Mx- in the low mass (low k±) region, we expect a non-perturbative picture of the 
type shown in Fig. 5 to give the dominant contribution. For this one-meson-exchange 
mechanism (see for instance [12,64,65]), the MM pair is created via double-Pomeron 
exchange, with an intermediate t-channel off-shell meson, and the amplitude is calculated 
using the tools of Regge theory, see [11] for more details 17 . In general, we should include 
both these contributions to the total meson pair CEP cross section (recalling from [11] that 
there is only a small interference between the dominantly imaginary perturbative and the 
dominantly real non-perturbative amplitudes), although in the limit of low/high meson k± 
only the perturbative/non-perturbative descriptions will be applicable. 

17 This subject has been also addressed in [12], although we have strong reservations concerning their 
treatment of the final state interaction between the outgoing mesons, for the reasons described above 
(5.14). 
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To calculate the non-perturbative MM contribution, we note that the CEP cross 
section is given by 

r CEP _ S 2 f j_ 2 jJ2 . j„. j_. ji.2 l-^P 



a = I67r(167r 2 ) 2 J d Pi±- d P2±dyidy 2 dk ± — , (5.11) 

where \fs is the c.m.s. energy, pi±,P2± are transverse momenta of the outgoing protons, k± 
is the meson transverse momentum and yi are the meson rapidities. S 2 is the soft survival 
factor - see below for a discussion of this. The matrix element is given by Ai = Mf + M«, 
with t = (Pi — k^) 2 , u = (Pi — k^) 2 , where Pj is the momentum transfer through Pomeron 
i, and £13,4 are the meson momenta. We have 

Mt = j^fM^M±)F 2 M(i)a 2 i^j (^J , (5.12) 

where M is the meson mass and we take sq = 1 GeV 2 and ap(p 2 i ) = 1.08 — 0.25|p 2 ± |, 
for p 2 ± measured in GeV 2 [98], and Sij = (p[ + kj) 2 is the c.m.s. energy squared of the 
final state proton-meson system (ij). The proton form factors are as usual taken to have 
an exponential form, F p (ti) = exp(Pj£j/2), while the slope of the Pomeron trajectory is 
included in the definition of the slope 18 , 



B % = 6 + 2a' log i^-j j (5.13) 

with bo = 4 GeV -2 . Concentrating on the case of tttt production, the overall cross section 
normalisation is set by the total pion-proton cross section a(irp) = ao(sij / so) a ^ _1 ~ 30 
mb at the relevant sub-energy for tttt production at the LHC. The factor F\j (t) in (5.11) 
is the form factor of the intermediate off-shell meson and, as discussed in [11], it is quite 
poorly known, in particular for larger values of t. It seems reasonable to take a typical 
'soft' exponential form Pm(*) = exp (b s (i — M 2 )), and the value of the slope can be 
approximately fitted to reproduce the correct normalisation of CERN-ISR data [99], as 
shown in Fig. 6, where a fair fit 19 to the data is given by the choice b Q g = 0.5 GeV -2 . 
While this value, which we note corresponds to quite a high average meson k\ ~ —t for 
a typical 'soft' process, reproduces the overall data normalisation, we note that the phase 
space region (M W7r < 2 GeV, \y w \ < 1.5) in which the data are collected places a very 
limited constraint on the cross section contribution with relatively large pion and so 
the higher t behaviour of Fm (t) , and therefore the non-perturbative cross section (which is 
sensitive to ~ F^), is still quite uncertain. With this is in mind we also show the prediction 
using the lower value of 6 flf = 0.625 GeV -2 , which we will also use when making numerical 
predictions in Section 6.1. While this seems to give too small a cross section compared to 
the CERN-ISR data, nonetheless we consider this value to indicate the sort of spread in 
predictions we may expect at higher pion k±, where the existing data place little constraint 



18 We note that in the equivalent formula to (5.12) in [11], the factors of a(pf ± ) should read ap(Q). 

19 The contribution of secondary Reggeons is also included here, using the fit of [100]. We are grateful to 



Piotr Lebiedowicz for providing a scan of the data found in [99] . 
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on the form of F n . 

It is worth emphasising that the theory curves shown in Fig. 6 should be considered only 
as a low (non-resonant) limit of the actual experimental cross section. The region < 2 
GeV manifests a number of resonances (/ 2 (1270), / (1370), / (1500), /£(1525), / 2 (1950), ...) 
which overlap with each other. Besides this, the t-channel state M* could correspond 
not only to pion exchange but also to the exchange of heavier states (a±, 02,...), which 
may enlarge the CEP cross section and explain the rather low value of the slope b Q s 
used to describe the data in Fig. 6 (see also [12]). Recall also that within the quasi- 
eikonal formalism the probability of higher mass pion excitation is about 0.45 [85], which 
is rather large, compared to 0.15 for the proton dissociation case, indicating a non-negligible 
contribution from high mass meson, M*, states 20 . 

In Fig. 7 we show the prediction using the same non-perturbative model for K + K~ 
CEP, compared to ISR data [102]. Far enough above the KK threshold, we expect the 
K + K~ differential cross section da/dMxK to be roughly suppressed by a factor 21 of 
(c/fp/oVp) 4 ~ 1/2 (at ISR energies this is complicated by the non-negligible contribu- 
tion of secondary Reggeons, although the level of suppression is still roughly the same). 
There will be some additional suppression from the Kaon mass tuk > m n nearer the 
KK threshold, and in particular the integrated K + K~ cross section will be signficantly 
suppressed due to this threshold effect. 

Finally, we have to include an additional suppression factor to calculate the genuinely 
exclusive MM (= ittv,...) cross section, i.e. that due to screening corrections. This was dis- 
cussed in [11], and here we give some clarifying remarks. We are considering the screening 
corrections to the non-perturbative amplitude shown in Fig. 5; in terms of the Reggeon 
formalism these corrections are described by the exchange of additional (one or more) 
Pomerons. First, there is the exchange between the two incoming (outgoing) protons 
(pi,P2), which is just the usual eikonal survival factor S^. Next, we have to consider the 
exchange between the upper (lower) proton and the lower (upper) meson, that is between 
Pi (P2) and M4 (Ms). We do not consider the exchange between p\ (p 2 ) and M3 (M4) 

20 On the other hand the effect of low mass proton dissociation, which is included in the models of [73,74] 
tends to reduce the survival probability, S^, an d therefore the exclusive cross section. In Fig. 6 we use 
these models to calculate the 'eikonal' survival factor, S%^. While there is tension between these models 
and the latest TOTEM data on elastic pp scattering at 7 TeV [101], with the models in particular giving a 
probability of low mass proton dissociation which is smaller than that indicated by the data, the updated 
model discussed in Section 6.1, which accounts for this, gives a cross section normalisation that is too low 
to fit the ISR data for realistic values of bog. We therefore treat the results of the one-pion-exchange model 
with this higher value of Se ik as giving a phenomenological description to the data, but note that a more 
complete treatment, including for example the exchange of other higher mass resonances (01,02,...) may 
be required to fully describe the ISR data. It is also not completely excluded that the higher value of «Sf ik 
that the previous models of [73,74] give may be more suitable, see Section 6.1 for further discussion. 

21 This observation is in strong contrast to the much higher level (~ 2 orders of magnitude at M„( K ^) = 4 
GeV) of suppression predicted in [103], see e.g. Fig. 7 (Right, dashed line), compared to Fig. 5 (Right, 
solid line) of [100] (Figure numbers refer to published versions) , which to our understanding are calculated 
using a very similar model to the one we describe here, with these plots in particular corresponding to the 
situation prior to the inclusion of -ktt(KK) rescattering, which may in principle complicate the issue. We 
can find no reason for a such an extreme level of suppression, and indeed Fig. 7 of [103] appears to be 
in tension with the ISR data [102]. Although contact was made with the authors to clarify this, no clear 
justification for this (very large) discrepency was given by them. 
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Figure 6: Differential cross section dcr/dM^ for the non-perturbative contribution to ir + ir~ CEP 
at *Js = 62 GeV, compared to CERN-ISR data [99]. The theory curves are calculated using the 
model described in the text, and with secondary Reggeons included, using the fit of [100]. In all 
cases, the pions are restricted to lie in the rapidity region \y^\ < 1.5 and the cut > 0.9 is 
imposed on the outgoing protons. 
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Figure 7: Differential cross section da/dMKK for the non-perturbative contribution to K + K~ 
CEP at y/s — 62 GeV, compared to CERN-ISR data [102]. The theory curves are calculated using 
the model described in the text, and with secondary Reggeons included, using the fit of [103]. In 
all cases, the kaons are restricted to lie in the rapidity region \yx\ < 1-5 and the cut > 0.9 is 
imposed on the outgoing protons. 



as for the corresponding pion-proton amplitude we have used the Donnachie-Landshoff 
parametrisation [98], in which the effective Pomeron Pi {P2) exchange already includes the 
sum of eikonal-like multi-(bare) Pomeron diagrams. Together with the possible exchange 
between the upper (lower) proton and the lower (upper) Pomeron, this corresponds to the 
enhanced screening, 5 en h, shown in Fig. 5. Note however that due to the rather small phe- 
nomenological value of the triple-Pomeron coupling (see for instance [86,91,92]) the main 
effect is expected to be from the secondary proton-meson interaction. However, in this 
case at higher values of V§ the size of the screening effect will be strongly suppressed by 
the small size of the produced pions ~ 1 /y/§. In the perturbative approach, this screening 
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effect is therefore power suppressed by the small absorptive cross section ~ 1/&j_ °^ the 
meson in the early stage of its formation 22 , and can be safely ignored in the phase space 
region we will consider here. In the nonperturbative case, unfortunately, there is no clear 
prescription for treating this rescattering of secondaries, as the cross section is written 
in terms of 'effective' Pomeron trajectories which are already renormalized by absorptive 
effects. As we do not know the details of how this renormalization occurs, it is therefore 
difficult to evaluate the level of proton-meson rescattering that we will expect to occur. 
Nonetheless, for reasonable values of the meson k± which we will be considering we expect 
the size of this rescattering to be similarly suppressed in the non-perturbative case, and 
we will ignore such an effect in our numerical calculations. 

Next, we recall that there is no secondary interaction between two outgoing mesons 
mediated by Reggeon exchange. This can be seen as follows: the time, r, needed to 
form a Reggeon increases linearly with energy, r oc E, and at high meson-meson energy, 
E ~ y/ s MMi ^ e t wo mesons travel apart at very close to the velocity of light. Since the 
meson pair production time in Fig. 5 is practically instantaneous (~ 1/E), while a much 
longer time (~ E/m 2 ) is need for the formation of a Reggeon by the secondary meson, 
there is insufficient time for a Reggeon emitted by one meson to interact with the other 23 . 
In other words, the two outcoming high energy mesons cannot talk to each other. Formally, 
when we consider the diagram with Reggeon exchange between the two outgoing mesons, 
with the Reggeon represented by a ladder-type graph, and check the imaginary part, ie, of 
the propagators, we see that the two (or three) poles which give the leading contribution at 
large s M j^ are all on the same side of the integration contour, and so closing the contour on 
the other side gives a zero leading s m jj contribution. As was shown by Mandelstam [107], 
the leading s contribution in the case of additional Reggeon exchange only comes from 
non-planar diagrams and not from planar graphs, of the type discussed in [108]. Thus the 
only possibility to have rescattering of the objects at high s m jj without violating causality 
is to consider the rescattering of incoming pomerons Pi and P2, but this is suppressed 
numerically by the small value of the triple-Pomeron vertex (see for instance [86,91,92]). 
We do not account for such an effect here. Instead, following [11], we introduce an extra 
suppression of the form of exp(— n), corresponding to the small Poisson probability not to 
emit other secondaries in the PIP — > M3M4 process at the initial meson pair production 
stage (rather than being due to final-state interactions between the mesons). Here n(s M jj) 
is the mean number of secondaries. This factor may be described as the Reggeization of 
the M* meson exchange, which means that we now deal with non-local meson-Pomeron 
vertices and the t-channel meson M* becomes a non-local object, i.e. it has its own size. 
It is this non-locality that is responsible for the non-violation of causality. More precisely 



22 To understand this we should recall that in the early production stage (in particular, in the case of a 
large f-channel meson M* virtuality) the meson colour field does not have sufficient time to regenerate, and 
the meson exists as a half-dressed (see e.g. [f04, 105]) particle (or pre-hadron, see [f06] for a recent review) 
with transverse size (much) smaller than that of the fully formed hadrons. Therefore, its interaction cross 
section should be weaker than in the case of an asymptotic (completely hadronized) meson state. 

23 We thank Yuri Dokshitzer for discussions about this space-time argument. 
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for the case of tttt CEP we take 



n(s 7r7r ) = 



M n7T < M /2 (i270) 




> M /2(1270) , 



(5.14) 



with sq = M\ 



/ 2 (1270) 



and c = 0.7- these precise values are taken for definiteness, but we note 



that in either case different choices are certainly possible. As discussed above, this choice 
appears to give a fair description of the existing ISR data. In this way, we account for the 
fact that we expect no additional suppression in the lower mass resonance region, where 
additional interactions at the meson pair production state lead mainly to the formation 
of resonances and not to the production of new secondaries. A similar although slightly 
modified procedure is taken for the case of K + K~ CEP. In particular we replace Mj 2 ( 12 7o) 
in (5.14) with Mj 2 ( 1525 ), and account for the fact that we should expect the number of 
secondaries h,(skk) to be a function of the free energy, Ef ree = M^k — 2M^, available 
for the creation of secondary particles, which can be numerically important because of the 
larger Kaon mass Mk- 

6. Results: predictions for the LHC, RHIC and the Tevatron. 
6.1 Meson pair CEP 

In this Section we present some numerical predictions for meson pair CEP at the LHC 
and RHIC. For the sake of brevity, we only consider the case of tttt CEP, although similar 
prediction for, say, KK CEP can readily be made (see Section 5.2). For the perturba- 
tive contribution, i.e. that found using the formalism outlined in Section 2, we use the 
MSTW08L0 PDFs, as these give a 77 CEP cross section that is in better agreement with 
the recent CDF data (see Section 3). For the non-perturbative contribution, we use the 
formalism described in Section 5.2. To calculate the soft survival factor we use the models 
of [73,74], incorporating new results for S 2 ik using the most recent fit [109], which includes 
the TOTEM data on elastic pp scattering at 7 TeV [101]. As the previous models of soft 
diffraction [73, 74] gave a somewhat lower value for cr tot at 7 TeV than the value measured 
by TOTEM, this most recent fit requires an increased optical density 0(6*) of the proton, in 
particular at low impact parameter, b t , see [109] for more details. Because of this increased 
optical density the probability of additional soft interactions is increased and the survival 
factor S^ ik calculated using this model is lower, giving a factor of ~ 2 decrease in the the 
combined S 2 = S^S 2 ^ for 77 CEP at y/s = 7 TeV, when compared with the previous 
predictions in [8-11,75]. However, at this stage it is possible that with further theoretical 
and/or experimental developments (for example, early measurements of diffractive dijet 
production by CMS [110] indicate that the gap survival probability may be higher than 
that given by our updated model), this value of S 2 ik may increase, and so the following 77 
CEP predictions should be viewed as slightly conservative with respect to the gap survival 
probability. 
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Figure 8: Perturbative and non-perturbative contributions to non-resonant 7r + 7r~ CEP differential 
cross sections dcr/dfc^ v at ^/s = 7 TeV, where fc^ v = (fc_i_(7r + ) + kj_(n~))/2 is the scalar average of 
the pion's transverse momentum, for the acceptance region k± > 0.5 GeV and \rj„\ < 1.0. 





(^/miii) '(max) 


A: ± (min)(|A:|(min)) [GeV] 


^Non-pert. [ nb ] 


ALICE 


(-0.9,0.9) 





840-1100 


ATLAS 


(-2.5,2.5) 


0.5 


1100-1700 


CMS 


(-2,2) 


0.2 


5300-6700 


LHCb 


(2,5) 


3 


3200-4000 


STAR 


(-1.0,1.0) 


0.15 


2500-3100 



Table 2: Pseudorapidity (rj) ranges and minimum k± values corresponding to the regions with 
good detection efficiency for tt + it~ CEP for the LHC detectors and the STAR detector at RHIC. In 
the case of the ALICE detector, an additional cut, Af W7r > 0.8 GeV, is imposed, and in the case of 
LHCb the momentum cut is placed on |/c|(min) and not k±. Also shown is the (non-perturbative) 
non-resonant tt + it~ CEP cross section, in nb, within these acceptance regions. The upper (lower) 
values of the cross sections correspond to taking b Q g = 0.5(0.625) GeV~ 2 , as described in Section 5.2. 



On the other hand, we find that for example the tttt CEP cross section is almost 
unchanged by taking this updated model for S?L. We recall that the new optical density 
£l(bt), motivated by the TOTEM data [101], increases mainly at low bt ■ On the other hand 
the LO gg —> tttt amplitude vanishes for J z = gluons, and so tttt CEP only occurs at LO 
for | J z \ = 2 fusing gluons. This contribution, due to its spin structure, h^btu, in impact 
parameter space, vanishes as bt — > 0. Therefore the more peripheral tttt CEP cross section 
is less sensitive (in comparison with 77 CEP where the J z = contribution dominates) to 
the increase in Q(bt) at low bt'- this is equally true for other CEP process (xc(i,2)> K + K~ ...) 
that vanish in the forward proton limit. 

In Table 2 we show the cuts that we will impose to roughly account for the pseudora- 
pidity and momentum acceptance of the LHC detectors and the STAR detector at RHIC 
for meson pair CEP [111]. We recall that the STAR detector is currently expected to have 
the ability to tag protons and therefore select purely exclusive events. At the LHC, with- 
out forward proton detectors exclusivity must instead be selected by vetoing on additional 
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k± > 1.5 GeV 


k± > 2 GeV 




Non-pert. 


Pert. 


Non-pert. 


Pert. 


ALICE 


23-120 


14 


0.058-0.69 


0.82 


ATLAS 


170-720 


50 


0.76-7.0 


3.0 


CMS 


130-580 


47 


0.53-5.2 


2.9 


LHCb 


69-310 


10 


0.23-2.2 


0.53 


STAR 


39-190 


1.1 


0.11-1.2 


0.067 


STAR ( P±i < 0.4 GeV) 


13-76 


0.31 


0.025-0.39 


0.018 


STAR {p ±i > 0.5 GeV) 


6.1-23 


0.15 


0.026-0.21 


0.008 



Table 3: Pcrturbative and non-perturbative contributions to non-resonant ir + ir~ CEP cross sec- 
tions (in pb) at the LHC (y/s = 7 TeV), for acceptance regions denned in Table 2, but with an 
additional cut on the pion k±, as indicated in the table. Also shown are the cross sections at 
RHIC energies (y/s = 500 GeV), with and without cuts on the outgoing proton transverse mo- 
menta p±i. The upper (lower) values of the non-perturbative cross sections correspond to taking 
b s = 0.5(0.625) GcV~ 2 , as described in Section 5.2 and in the text. 

hadronic activity in a large enough rapidity region, although in this case there will also 
be some contribution from non-exclusive events, in particular where one or both of the 
protons dissociate, as in (1.2). Such a background contribution could be particularly im- 
portant at LHCb, where the veto region (2 < r\ < 5 and — 4 < r\ < — 2 for charged tracks) 
is fairly limited, see Section 4 for further discussion. However, as discussed in the Introduc- 
tion, recently the possibility of measuring exclusive low mass pion pair production using 
the ALFA RP detectors during the special low pile-up LHC runs has been proposed [31], 
and data taken as part of a proposed dedicated CMS run of single no-pile-up interactions, 
engaged and taking data in common with TOTEM would be especially promising [48]. 

Using the cuts of Table 2 and taking 7r + 7r~ as a specific example 24 , we then also show 
in Table 2 the non-resonant tt + tt~ CEP cross section. As discussed in Section 5.2, the size 
of the non-perturbative cross section at larger values of the meson k± is largely uncertain, 
due to our lack of knowledge of the form factor of the intermediate off-shell pion F n . 
We therefore show approximate upper (lower) bounds on the range of predictions, found 
using the values for the slope parameter b Q g = 0.5(0.625) GeV -2 , assuming an exponential 
functional form F w (t) = exp (b g(t — M 2 )). In this case the predicted cross section, which is 
primarily sensitive to lower values of the meson k± , does not change too much between these 
possibilities, although as we shall see this is no longer true when k± cuts are imposed. We 
only show the cross-section prediction calculated within the non-perturbative framework, 
because at the low M nn and k±(ir) values accessible due to the reasonably low k± reach of 
the current LHC experiments and the relatively large rapidity coverage, the application of 
the pQCD approach used to describe both the overall CEP process, as well as the gg — > 
7r + 7r~ subprocess cross section, cannot be justified. In particular, it gives an unphysical 
divergent answer as the scale of as and the off-shellness of the propagators in the gg — > 
7r + 7r~ subprocess decrease. 

24 We recall that the iy°-k° CEP cross-section prediction is 1/2 of the iv + it~ cross section, from isotopic 
invariance and the identity of the final state w° states. 
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We can see in Table 2 that the expected n + ir~ cross sections, calculated within the 
non-perturbative framework, are quite large. This can even be the case as we increase 
Mx, when the pions are produced at non-zero pseudorapidity and can still have fairly- 
low k±. However, in Fig. 8, we show the differential cross section da/dkY (where k± 
is a suitably defined transverse momentum variable, see the Figure caption) in both the 
perturbative and non-perturbative case, for the two choices of b $ described above. Clearly 
there is a strong damping in the non-perturbative cross section at higher k±, caused by 
the soft pion form factor F n (k^_), which drops sharply with the pion k±: in this region, the 
high k± tail of the it + it~ CEP process is generated by a purely pQCD mechanism, which 
cannot be predicted within the framework of Regge theory. However we can see that the 
different choices of b a e, while predicting fairly similar total cross sections (see Table 2), 
give roughly an order of magnitude variation in the differential cross sections at higher k± . 
Furthermore, as discussed in Section 5.2, the existing CERN-ISR data [99] imposes some 
constraint on the off-shell pion form factor, F n , at low but not high pion k±, and so the 
non-perturbative cross section may in principle lie outside this range (although for such a 
'soft' process we would not expect the slope parameter to be much smaller than the lower 
choice b Q g = 0.5 GeV -2 ); these estimates represent our best educated extrapolation to the 
higher k± region. 

Nonetheless, if we impose suitable cuts on the pion k±, the perturbative contribution 
should dominate and can be used to reliably estimate the full 7r + 7r~ CEP cross section. In 
Table 3 we show cross-section predictions for the same experimental pseudorapidity regions, 
but with the additional cut on the pion k± > 2.0 GeV, and the perturbative mechanism is 
expected to be non-negligible and possibly dominant. For higher values of the k±_ cut (see 
Fig. 8) we therefore expect the non-perturbative contribution to be negligible, although 
we can see in Table 2 that this is in a region where the perturbative cross section (~ 1 
pb or lower) is quite small and unfortunately may not be experimentally accessible for the 
low luminosity runs suitable for CEP observations. This is due to the strong (~ 1/100) 
suppression by the J z = selection rule of the perturbative 7r + 7r~ CEP cross section 25 that 
comes from the fact that the gg — > 7r + 7r~ amplitude vanishes at LO for J z = incoming 
gluons [11]. We also show in Table 3 the irir CEP cross section with the cut k± > 1.5 GeV: in 
this case, while the cross sections are higher, the non-perturbative contribution is expected 
to be dominant. As noted in [11], other related processes, such as rj(')rj(') CEP, where the 
production amplitudes do not vanish for J z = incoming gluons, are expected to have 
significantly larger cross sections and so may represent experimentally more realistic and 
theoretically cleaner (that is, with a smaller non-perturbative contribution) observables. 

We recall that at the STAR detector at RHIC, a large CEP data sample with tagged 
forward protons is expected to be available in the near future [23-25]. The position of 
the existing and proposed proton tagging detectors is ideally suited for the observation 
of meson pair CEP in the experimentally accessible mass region. In Table 3 we also 

25 Recalling the observed enhancement by LHCb [27] of the Xc2 rate above theory expectations, which 
may be due to a significant contamination from non-exclusive events, we note that the observed tt + it~ rate 
may be correspondingly enhanced by such contamination at the LHC, where the outgoing protons are not 
tagged. 
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da/dM„ [pb/GeV], ^ft=7 TeV, \rf \ < 1 dc/AMm, [pb/GeV], v's = 7 TeV, |jf| < 1, k\ > 1.5 GeV 




Figure 9: \cQ — > 7r + 7r~ Breit-Wigner peak, and perturbative and non-perturbative (for different 
b s values) contributions to non-resonant tt + tt~ CEP in the XcO mass region at y/s = 7 TeV, for 
different and 7/ cuts on the pions. 



da/dM„ [pb/GeV], <fi = 1.96 TeV, < 1 da/dM„ [pb/GeV], yfi = 1.96 TeV, |jf| < 1, fcf > 1.5 GeV 




Figure 10: \c0 —> tt + tc~ Breit-Wigner peak, and perturbative and non-perturbative (for different 
b ff values) contributions to non-resonant ir + ir~ CEP in the XcO mass region at y/s = 1.96 TeV, for 
different k± and rj cuts on the pions. 



show the 7r + 7r~ cross section after acceptance cuts on the proton transverse momenta p±i 
have been imposed, with in particular the first cut {p±i < 0.4 GeV) corresponding to the 
RPs currently installed (Phase I) and the second cut (p±i > 0.5 GeV) to the phase II 
RPs, to be installed in ~ 2013. Such an evaluation requires an explicit inclusion of the 
proton distributions after the inclusion of soft survival effects. We can see that the 
relative perturbative to non-perturbative contribution is expect to be much lower than at 
\fs = 7 TeV, which is due to the different y/s dependence of the production mechanisms: 
while the perturbative contribution drops quite steeply with y/s due to the fairly steep 
x-dependence of the gluon density, g(x,Q 2 ), the non-perturbative contribution exhibits a 
less steep, Regge-like energy dependence, see (5.12), and may even increase with decreasing 
y/s due to the larger survival factor S^ lk . It therefore appears that at the STAR detector 
the perturbative contribution will be negligible, although we note that this mechanism, 
which will dominantly occur for \J Z \ = 2 pions, will exhibit a different distribution in <p, 
the azimuthal angle between the outgoing protons [8], than in the non-perturbative case, 
where no such restriction on the pion angular momentum is present. The observation of 
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d<j/dM„ [pb/GeV], = 0.5 TeV, \if\ < 1 da/AM m [pb/GeV], y/i = 0.5 TeV, \rf\ < 1, fcj > 1.5 GeV 




Figure 11: % c o — > -k + h~ Brcit-Wigner peak, and perturbative and non-perturbative (for different 
b a g values) contributions to non-resonant ir + ir~ CEP in the \c0 mass region at y/s — 0.5 TeV, for 
different k± and r\ cuts on the pions. 



7r + 7r CEP with tagged forward protons would therefore represent an interesting test of 
these two mechanisms and their relative importance. 





k± > 2 GeV 


k ± > 2.5 GeV 


k± > 5.5 GeV 


y/s = 1.96 TeV (\rj\ < 1) 


9.0 


3.1 


0.059 


v /i= 7 TeV (\rj\ < 2.4) 


19 


6.1 


0.12 



Table 4: Perturbative cross section ratio cr pcrt ' (ir + Tr~)/a(jj) (in %) at different c.m.s. energies 
and for different cuts on the pion/photon k± and rj. MSTW08LO PDFs are used, with the survival 
factors Sg ik , S^ nh calculated as described in the text. 



In Table 4 we show the (perturbative) cross section ratio a peTt ' (tt + tt~) / a(^^) for dif- 
ferent cuts on the pion/photon fcj_(7r, 7) at the LHC (y/s = 7 TeV): for such an observable, 
the theoretical uncertainties from the choice of PDF (and to a lesser extent those due to 
the survival factors S^ nh , S^ ik ) practically cancel out, and such ratios therefore provide a 
better test of the underlying theory. As we can see, for k± > 2 GeV the tt + tt~ cross section 
is predicted to be a factor of ~ 5 smaller than the 77 cross sections. Any excess above this 
might be evidence for a non-negligible non-perturbative 7T7t contribution, which Table 3 
indicates may be present; for the higher k± cuts in Table 4 such a contribution can on the 
other hand be safely ignored, see Fig. 8. As k±_ increases the 7r + 7r - cross section is strongly 
suppressed, which is due to the factor of {fir/k±) A present in the gg — > tt + tt~ cross section, 
representing the small probability for each outgoing qq pair to form a pion. We recall that 
the cross section ratios shown in Table 4 are calculated using the new results for from 
the most recent fit [109], which includes the TOTEM data on elastic pp scattering at 7 
TeV [101], and results in a factor of ~ 2 decrease in the 77 cross section, while leaving the 
more peripheral 7T7t cross section roughly unchanged. However, as described in the begin- 
ning of this Section, it is possible that this lower value of S^ ik is somewhat conservative, 
and so the cross section ratios shown in Table 4 may be slightly lower. For comparison, we 
also show the same cross section ratio cr pert '(7r + 7r~)/o"(77) at the Tevatron, but using the 
previous (higher) values of S 2 calculated using the models described in [73,74], and used 
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in for example [9,11], see Section 6.2 for further discussion. 

Finally, in Figs. 9-11 we show results for the tt + it~ CEP cross section in the XcO 
mass region at LHC, Tevatron and RHIC energies, which will be relevant for evaluating 
the potential continuum background to resonant XcO — > tt + tt~ production. The it + it~ 
mass distribution from XcO decay is given by a simple non-relativistic Breit-Wigner, with 
the XcO cross section normalisation set using the Super CHIC MC [58], which implements 
the theory described in [8]. We can see that once basic r] cuts are imposed on the final 
state pions, the XcO signal is expected to lie at a similar level to the non-perturbative 
background with b R = 0.5 GeV -2 , and that the perturbative contribution is expected to 
be negligible. Although we note that the non-perturbative background may be somewhat 
lower in this mass region (for comparison we also show the background for the choice 
6 off = 0.625 GeV -2 , which gives a lower cross section), it is not completely clear that the 
signal peak will be visible over the background, taking into account the various theory 
uncertainties (experimental resolution effects may also decrease the S/B ratio). However, 
if we also impose a simple k± > 1.5 GeV cut on the final-state pions, we can see that 
the background is strongly reduced with a much smaller effect on the Xc signal rate (for 
which the Xc mass M x ~ 3.5 GeV ensures that a large fraction of the central pions have 
k± > 1.5 GeV): the predicted XcO tt + tt~ rate lies (at least) an order of magnitude 
above the expected background. We can therefore safely conclude that even within the (in 
principle quite large) theory uncertainties, XcO — > vr + 7r~ is expected to represent a clean 
experimental signal, with a low continuum background once suitable cuts are imposed. 

6.2 77 CEP: PDF comparison 





MSTW08LO 


CTEQ6L 


GJR08LO 


MRST99 


CT10 


NNPDF2.1 


s/s = 1.96 TeV (|r?| < 1) 


1.4 


2.2 


3.6 


0.35 


0.47 


0.29 


Vi = 7TeV (M < 1) 


0.061 


0.069 


0.16 


0.013 


0.0094 


0.0057 


^ = 7 TeV (\r]\ < 2.5) 


0.18 


0.20 


0.45 


0.039 


0.027 


0.017 



Table 5: 77 CEP cross sections (in pb) for different choices of gluon PDF, at yfs = 1.96 and 7 
TeV, and for different cuts on the photon pseudorapidity, rj. The photons are restricted to have 
transverse energy E± > 2.5 GeV at y/s = 1.96 TeV and E± > 5.5 GeV at y/s = 7 TeV. 



As discussed in Section 3.1, even after the various theoretical (and experimental) un- 
certainties are taken into account, there is some tension between the recent CDF 77 CEP 
data [19] and the theoretical prediction using the NLO MRST99 PDF set, while if the 
MSTW08LO PDF set is taken there is good agreement between data and theory. Al- 
though some caution is needed, recalling the theoretical uncertainties, we can in principle 
use these CEP measurements to shed some light on the gluon PDF in this low-x and low- 
Q 2 region, where it is poorly determined. We recall in particular the significant difference 
between the LO and NLO PDF fits: while the LO PDFs have quite a steep low-x depen- 
dence, the NLO PDFs are much smaller and can even be negative, at small x and Q 2 for 
the more modern fits. With this in mind, in Table 5 we show cross section predictions 
for 77 CEP within the CDF kinematics, at = 1.96 TeV, and we can see that for a 
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representative sample of LO PDF sets, there is good agreement with the data, while the 
NLO set predictions (we are careful to choose only those sets which are not negative in the 
relevant kinematic region) all lie significantly below it. Future LHC measurements would 
allow us to begin to study the energy dependence of the 77 cross section, and through this 
we may hope to provide further information about the gluon PDF, in particular at even 
lower values of x. In Table 5, we therefore also show predictions for the LHC (y/s = 7 
TeV) with the same PDF sets, and for different cuts on the photon pseudorapidity, n. 

We note that the values in the lowest row (\t]\ < 2.5) are somewhat lower than those 
given in Fig. 7 of [63], which is principally due to the different (lower) value for S 2 ik we 
now use to account for the recent TOTEM data [101], as well as various other theoretical 
refinements (the 'enhanced' survival factor, S 2 nh , non-zero proton pj_ etc) not included in 
the ExHuME MC [112] that is used to calculate these numbers 26 . On the other hand, as 
discussed at the beginning of Section 6.1, we note that the value of S 2 lk may in principle 
be higher for 77 CEP, in which case the LHC cross section predictions in Table 5 will also 
be correspondingly increased. With this in mind, for the Tevatron predictions we keep 
the previous (higher) values of S 2 calculated using the models described in [73,74], and 
used in for example [9, 11], consistently with the cross section predictions quoted in [19]. 
However, we note that using the value of S 2 ik (S 2 nh ) calculated using the fit of [109] the 
corresponding 77 cross section prediction will be a factor of ~ 2 lower, although this 
appears to be disfavoured by the CDF data. 

7. Conclusions 

We have performed a detailed theoretical and phenomenological study of several bench- 
mark central exclusive production processes in high-energy proton-(anti)proton collisions, 
focusing in particular on 77 and meson pair CEP at the Tevatron, RHIC and the LHC. 
Our approach is based on the successful and well-motivated Durham perturbative model, 
but we have also addressed the importance of additional non-perturbative contributions 
in the case of meson pair CEP, which are expected to be important for lower meson k±. 
Unlike the case of inclusive production of massive particles and final states, even in the 
perturbative sector there are large uncertainties in the predictions, in particular from un- 
known higher order perturbative corrections, from the choice of gluon PDFs and from the 
soft survival factors. 27 

Our current understanding of these uncertainties can be summarized as follows. 

(i) The dependence on the gluon PDF is amplified by the fact that the CEP cross section 
is essentially proportional to (xg(x)) A . We have argued that it is reasonable to take 
the range spanned by leading order and next-to-leading order PDFs as indicative 
of the uncertainty, and have used MRST08 LO and MRST99 NLO for this purpose. 

26 We should also recall that, as discussed in Section 4, accounting for the possibility of proton dissociation 
would lead to an increase of the theoretical CEP prediction by about a factor of 2, which is consistent with 
the CMS upper limit. 

27 In practice it is only approximately true that we can separate these effects, since the emission of 
additional gluons can contribute to all three. 
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We noted that the former (which has a higher density of small-x gluons) gives better 
agreement with the recent CDF diphoton data [19]. We emphasise again that at the 
low x values that are relevant for the cases of charmonium, diphoton and dimeson 
CEP, the difference between the LO and NLO PDFs can be very large. 

(ii) Our CEP 'hard' matrix elements are calculated at leading order in QCD perturbation 
theory. NLO (and higher order) corrections come from a number of sources, of which 
the most important are: 

(a) NLO corrections to the hard gg — > X matrix element. 

(b) NLO effects in the unintegrated parton densities, see [83]. 

(c) Effects of radiation off the screening gluon, of the type discussed in [84]. 

(d) A systematic account of self-energy insertions in the propagator of the screening 
gluon 28 , see the discussion in [80]. 

Interestingly, all of the NLO effects identified above indicate a possible increase in 
the rate predicted at LO, by a total factor of order ~ 2, although we should recall 
that such effects cannot be considered completely seperately from the (potentially 
large) variation in the cross section which can occur when different (LO and NLO) 
PDF sets are used. 

(iii) Soft survival factors (both eikonal, S" 2 ik , and S^ nh ). Within the Durham soft diffrac- 
tion models [50, 51, 73] their overall effect is fairly stable, fluctuating by roughly 
30 — 35% for different sets of soft model parameters. 

(iv) For meson pair CEP, we also need to consider the impact of non-perturbative contri- 
butions, for example where the meson pair is created by double Pomeron exchange. 
These are only likely to be important for low invariant mass and/or transverse mo- 
mentum final states, and can therefore be suppressed by imposing suitable cuts on 
the meson k±. Nonetheless, for example for tttt CEP, where the perturbative contri- 
bution is strongly suppressed by the J z = selection rule which operates for CEP, a 
cut of k± > 2 GeV (or possibly even higher) is expected to be needed to safely sup- 
press the predicted non-perturbative contribution. In this case, the corresponding 
perturbative cross sections are of order ~ 1 pb or lower, and so this may be difficult 
to access experimentally. Other observables such as vOvC) CEP are expected on the 
other hand to have much larger perturbative cross sections and therefore relatively 

28 Usually in Feynman diagrams we should include a wavefunction renormalization factor Z 1 * 12 for each 
external parton leg. However the conventional PDFs are defined in a different way, with the whole Z 
factor (for each of the two legs) included in the PDF, while the 'hard' matrix element is calculated by 
taking Z = 1 for each of the external parton legs. For standard inclusive processes this prescription 
leads to the correct result, but some care is needed in the CEP case. Since the soft screening gluon in 
Fig. 1 is not connected to the hard matrix element, the left gluon propagator between the unintegrated 
skewed distributions f g {x\, ...) and f g (x2, ...) will be assigned a factor of Z 2 , instead of the usual Z. The 
corresponding first order a s correction therefore corresponds in part to an NLO contribution, which should 
not be included at lowest order in a consistent treatment. 
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smaller non-perturbative contributions. We also note that the observation of K + K 
(pp, AA...) CEP would also be of interest. 

To summarize, if we take into account all of the main sources listed above, we estimate 
the uncertainty in the CEP cross section predictions (on top of the choice of PDF set) 
to be of order ~ *2 for the case of 77 CEP. We would expect a larger uncertainty in 
the 7T7T predictions (even after the imposition of cuts to suppress the non-perturbative 
contributions), due to the strong suppression in the irn CEP cross section from the J z = 
selection rule, as in this case higher order/twist effects which allow even quite a small 
contribution from a J z = component may enhance the rate. Cross section ratios, in 
particular (1(77) / 'o{jtit) (as well as o-(x c )/o'(TnT) and cr(xb)/ a (ll) m t ne relevant mass 
regions), should be much less sensitive to the choice of PDFs and the survival factors, S 2 , 
and should therefore provide a good test of the perturbative approach to CEP. 

We have also discussed the importance of proton dissociation contributions. The stan- 
dard CEP formalism assumes that the elastically scattered protons are detected in the final 
state. So far, however, the only experimental access to such 'exclusive' events has been via 
large rapidity gap triggers. There is then the issue of the importance of proton dissociation 
contributions to the cross section. If the experimental configuration restricts the mass of 
the dissociated state to a few GeV, for example by vetoing on additional particle produc- 
tion over a wide rapidity range, then the cross section enhancement factor is likely to be 
of only a few tens of %. The effect of high-mass dissociation is much harder to estimate, 
but such events may be suppressed by restricting the transverse momentum of the central 
state. 

We finish by recalling that all the attempts to evaluate CEP quantitatively (for suffi- 
ciently high central system mass Mx) are based on models which contains the same basic 
elements as those in the Durham model discussed in this paper. Here we consider the 
theoretical uncertainties contained in the various theory elements and the constraints that 
come from a phenomenological analysis of the existing CEP data. We have tried to make 
our analysis as full and complete as currently possible, although we hope in the future 
that further theoretical progress will be made, guided by new experimental data. It will in 
particular be very interesting to compare future RHIC and LHC data with our predictions 
for the processes discussed in this paper. For example, the forthcoming CMS data on 77 
CEP should help to reduce these uncertainties and should in particular clarify the question 
of the most appropriate choice of PDF set [63] . 
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